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^ ■ Abstract 

The effect of the on-ramp and off-ramp positions ii and ^2, respectively, on the one 

^ I dimensional-cehular automaton traffic flow behaviour, is investigated numericahy. The 

■ on-ramp and off-ramp rates at ii and 12 are ao and Pq, respectively. However, in the 
00 . open boundary conditions, with injecting and extracting rates a and /? and using par- 
^ I allel dynamics, several phases occur, depending on the position of ii by respect to 22- 

' Namely, low density phase (LDP), intermediate density phase (IDP), plateau current 
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+^ I ticritical behaviours take place in the (11,00) phase diagrams. 
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1 Introduction 



The field of transport have attracted increased attention in recent years [1-4]. This 
interest is due primarily to the fact that transportation problems are related to the 
global behaviour of systems with many elements interacting at short distances, such 
as the vehicles travelling on the streets or informations which travel over the internet 
network. 

In particular, the investigation of open traffic systems with on- and ofF-ramps is quite 
popular at the moment [5-10]. One reason for this is the impact of the understanding 
of varying the different flow rates in order to optimise the total flow or trip times. 
Among the different methods of investigation and simulation of highway traffic flow, 
asymmetric simple exclusion process (ASEP) is the most promising [11-13]. Indeed, 
ASEP is the simplest driven diffusive system where particles on a one- dimensional lat- 
tice hop with asymmetric rates under excluded volume constraints. 
The effect of on- and off-ramps have been largely studied by several researchers in the 
traffic flow field [5,11]. Our principal aim in this paper is to bring an explanation on the 
following point: let us suppose that we have a highway running in an urban conglomera- 
tion and that there are an access from urban conglomeration to the highway and an exit 
from the highway. We want to understand where the access and exit positions must be 
localised in order to maximize the flux of cars in the road. In other words, we want to 
study the effect of the on-ramp and off-ramp positions on the one dimensional-cellular 
automaton traffic fiow behaviour, in the open boundaries case. To this end, we have 
established the phase diagrams in the (ii, ao)-plane. Depending on the injecting and 
extracting rate values, an adequate localization of the on- and off-ramp positions leads 
to the appearance of new phases, topologies and several multicritical behaviours. More- 
over, to compare our results to those where only one off-ramp was taken into account 
[5], quantitative differences can be understood from the behaviour of average density 
and phase diagrams for different parameters. 

The paper is organised as follows: Model and method are given in section 2; section 
3 is reserved to results and discussion; the conclusion is presented in section 4. 
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2 Model and method 



We consider a one-dimensional lattice of length L. Each lattice site is either empty or 
occupied by one particle. Hence the state of the system is defined by a set of occupation 
numbers ri,r2,...,Tj^, while Tj = 1 (r^ = 0) means that the site i is occupied (empty). We 
suppose that the main road is single lane, an on-ramp and an off-ramp connect the main 
road only on single lattice ii for entry and on single lattice 12 for way out. During each 
time interval At, each particle jump to the empty adjacent site on its right and does not 
move otherwise (i 7^ ^2). At is an interesting parameter that enables the possibility to 
interpolate between the cases of fully parallel {At = 1) and random sequential {At — > 0) 
updates [12]. Particles are injected, by a rate aAt, in the first site being to the left 
side of the road if this site is empty, and particles enter in the road by site ii, with a 
probability agAt without constraint, if this site is empty. While, the particle being in 
the last site on the right can leave the road with a rate (5At and particles removed on 
the way out with a rate PoAt. At site ii {12) the occupation (absorption) priority is 
given to the particle which enter in the road (particle leaving the road). Hence the cars, 
which are added to the road, avoid any collision. 

If the system has the configuration Ti{t), T2{t),...,TL{t) at time t it will change at time 
t -I- At to the following: 
For i = ii , 

Ti{t + At/2) = 1 (1) 

with probability 

Qi = Ti{t) + K(l - ri{t)) - Ti{t){l - Ti+l{t))]At (2) 

and 

Ti{t + At/2) = (3) 

with probability 1 — qi. Where ii and cto denote the position of the entry site and the 
injection rate, respectively. 
For i = 12 , 

Ti{t + At/2) = 1 (4) 



with probability 
and 



Qi = Ti{t) + [Ti_l(l - n) - PoTi{t)]At (5) 



Ti{t + At/2) = (6) 



with probability I — Qi- Where 12 and /3o denote the position of the absorbing site and 
the absorbing rate, respectively. 
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For l<i<L with i ^ ii and i ^ 12, 

Ti{t + At) = l (7) 

with probabiHty 

Qi = nit) + - Ti(t)) - Ti{t){l - ri+i(i))]Ai (8) 

and 

Ti{t + M) = Q (9) 

with probability I — Qi- 
For 1 = 1, 

ri(i + At) = l (10) 

with probabiHty 

qi = Ti{t) + [a{l - n{t)) - n{t){l - T2mAt (11) 

and 

Ti{t + At) = (12) 

with probability 1 — gi- 
For i = L, 

TL{t + At) = l (13) 

with probability 

QL = Mt) + [TL-lim - rL{t)) - PTL{t)]At, (14) 

and 

TL{t + At) = (15) 

with probability 1 — qr,. 



3 Results and Discussion 



In our numerical calculations, the rule described above is updated in parallel, At = 1, 
i.e. during one update step the new particle position do not influence the rest and only 
the previous positions have to be taken into account [13]. During each of the time steps. 
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each particle moves one site unless the adjacent site on its right is occupied by another 

particle. 

In the following, the temporal average of any parameter is computed for 5 x 10^ to 10^ 
time steps. Starting the simulations from random configurations, the system reaches a 
stationary state after a sufficiently large number of time steps. In all our simulations, 
we averaged over 60 — 100 configurations. For the update step, we consider two sub 
steps as has been used in a previous work [14]. 

As we have mentioned previously, our aim in this paper is to study the effect of the 
on- and off-ramps positions ii and 12, respectively, for different values of a, (3q and 
on the average density and flux in the chain. The study is made in the open boundary 
conditions case. The length of the road studied in this paper is L=1000. Moreover, our 
goal here is to study the cases where ii is upstream or downstream 12, we then fixed 
12 = 500 and varied ii. 

The figure 1 gives the variation of the average density p as a function of the inject- 
ing rate qq for several values of the on-ramp position ii and for a = 0.1, /? = 0.1 
and /3o = 0.4. However, when the on-ramp is located upstream the off-ramp, the sys- 
tem studied exhibits four phases, namely: i) The low density phase (LDP), where the 
average density increases when increasing the rate of injected particles ao- ii) The in- 
termediate density phase {ID Pi) characterised by a smoothly increase of the density, 
iii) The plateau current phase {PC Pi) for which the density and current are constant 
in a special interval of ag- iv) The high density phase {HDPi) in which, for high values 
of ao, the current decreases and the density reaches a large value and remains constant. 
On the other hand, when ii is located downstream Z2, the high density phase disappears. 
Note that we add index 1 to each phase when this one is specific to the case ii < 12 and 
index 2 in the contrary case. In addition, when increasing ii for a given value of oq, the 
average density remains constant in the LDP, while it decreases in the IDPi and PCP 
then increases in the HDPi. Moreover, an inversion point is located at a special value 
of ao = aoi which corresponds to the PCPi — HDPi transition. This inversion can be 
explained as follows: 

- For ii < 12: if ao < aoi and ii increases, the section between ii and i2 decreases re- 
ducing the average density (a is small). While, if ag > aoi, there is a fast accumulation 
of particles in the section between ii and L {(3 being weak) which increases the average 
density. 

-For ii > 12: if ao is small, the increase in ii reduces the zone ii — L thus decreasing 
the average density. Whereas, when ao is large, the accumulation of particles between 
ii and L is compensated by the reduction in the section ii — L. The average density 
remains thus constant. 

In addition, if ii < 12, we note that the IDP, which doesn't appear in the model where 
only the off-ramp is taken into account [5], occurs for the intermediate values of ao 
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(aoci < ao < aoc2)- ooci and aoc2 correspond to the transition between LDP — IDP\ 
and IDPi — PC Pi, respectively. While, the PC Pi arises between two critical values 
aoc2 and aoc3 of injecting rate uq, where aoc3 corresponds to the PC Pi — HDPi transi- 
tion, which coincide with aoi- We note that this transition disappears when ii is located 
after 12 (Figure 1). 

In order to have a suitable criterion for determination of the nature of the transition, we 
identify the first order transition (abrupt transition) by a jump in the average density 
or by the existence of a peak in the derivative of p(ao) with respect to ao [13]. This 
means that the transitions shown in figure 1 are of first order type. 
Collecting the results illustrated in figure 1, the different regions are given on the phase 
diagram (ii,ao) shown in Figure 2. Beside this, such phase diagram exhibits the critical 
and multicritial behaviours. The critical end-points, around which there is no distinction 
between the phases, are denoted by CEP. In order to describe the different entities in 
the phase diagrams, the Griffiths notations [15] will be adopted: the multicritical point 
j\rnj^n jgnotcs the intersection of m lines of first order and n lines of second order. 
The indicates then the triple point, which is the intersection of three lines of first 
order. The critical point, which is the intersection of a line of first order and a line 
of second order, is denoted by C. However, the figure 2 shows two multicritical points 
A^B and five single critical end-points CEP. Indeed, for a given value of ao, the tran- 
sitions ID Pi — PCP2 and PC Pi — IDP2 are of the second order, while the transition 
HDPi — IDP2 is of the first order. We identify the second order transition by respect 
to the average density, which is continuous but presents an angular point [12]. For low 
values of a, [5 and /3o {a = 0.1, fi = 0.1 and /?o = 0.1), the (ii,ao) phase diagram is 
presented in figure 3, where all the transitions are of the first order, except that located 
between PCP1—IDP2. Moreover, such figure exhibits three critical end-points and one 
multicritical point of type A^B. The comparison of figures 2 and 3 shows up the effect 
of /3o on the (ii, ao) phase diagram. Indeed, for intermediate value of Pq, the IDPi and 
PCP2 take place. 

Furthermore, for a sufficiently large value of Pq {(3q = 0.8, as shown in figure 4), the 
IDPi disappears, two critical points (C), two triple points {A?) and two end-points 
(CEP) occur. 

Finally, to highlight the effect of /3, we present on figure 5 the (ii, ao) phase diagram for 
a = 0.1, /3o = 0.1 and (3 = 0.3. In this case, the system exhibits one multicritical point 

A^B, one critical point C, one triple point and two critical end-points CEP. From 
figures 2, 3 and 5, we deduce that the IDPi arises for intermediate values of extracting 
rates /3 or /3o, in agreement with the previous work [14]. 
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4 Conclusion 



Using numerical simulations, we have studied the effect of the on- and off-ramp positions 
on the traffic flow behaviour of a one dimensional-cellular automaton, with parallel up- 
date. When ii is upstream ^2, an inversion point takes place in p{ao)j PC Pi — HDPi 
transition. Depending on the values of a, [5 and /3o, the (zi,Q!o) phase diagram exhibits 
different topologies. The transitions between different phases are of the first or second 
order. Furthermore, the system exhibits a multicritical, critical and critical end-points. 
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Figure captions: Fig. 1: Average density p versus the injection rate uq for a = 0.1, 
/? = 0.1, /?o = 0.4. The number accompanying each curve denotes the value of ii. Fig. 
2: Phase diagram (ii,ci;o) for a = 0.1, /3 = 0.1, Po = 0.4. The solid and dot lines 
indicate the first and second order transitions, respectively. 
Fig. 3: Phase diagram (ii,ao) for a = 0.1, Pq = 0.1 and (3 = 0.1. 
Fig. 4: Phase diagram (ii,Q!o) for a = 0.1, Pq = 0.8 and P = 0.1. 
Fig. 5: Phase diagram (ii,ao) for a = 0.1, Pq = 0.1 and P = 0.3. 



8 



1,0 




1000 



800- 



PCP, 



600- 



a'b 



CO 

o 

Q_ 
Q_ 

E 



O 



400- 



□_ 
Q 



/ 

CEP 



200- 



IDP 







0,0 



0,2 



0,4 



IDP, 



0,6 



a'b 





HDP^ 


PCP. 





CEP 



0,8 



1,0 



a. 



1000 



800- 



IDP. 



CO 

o 

Q_ 
Q_ 

E 



600- 



400- 



CL 
LU 
O 



a'b 



o 



200- 



LDP 



PCP 



HDP 







0,0 



0,2 



0,4 



0,6 



0,8 



1,0 



a. 



1000 



800- 



600- 



O 

Q_ 
Q_ 

E 



O 



400- 



200- 







0,0 



PGP, 



LDP 



0,2 



CEP 



IDP. 



CEP 



PCP 



HDP 



0,4 



0,6 



0,8 



1,0 



a. 



1000 



800- 



600- 




400- 



200- 



